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We examine the relaxation of the Kob- Andersen Lennard- Jones binary mixture using Brown- 
ian dynamics computer simulations. We find that in accordance with mode-coupling theory the 
self-diffusion coefficient and the relaxation time show power-law dependence on temperature. How- 
ever, different mode-coupling temperatures and power laws can be obtained from the simulation data 
depending on the range of temperatures chosen for the power-law fits. The temperature that is com- 
monly reported as this system's mode-coupling transition temperature, in addition to being obtained 
from a power law fit, is a crossover temperature at which there is a change in the dynamics from 
the high temperature homogeneous, diffusive relaxation to a heterogeneous, hopping-like motion. 
The hopping-like motion is evident in the probability distributions of the logarithm of single-particle 
displacements: approaching the commonly reported mode-coupling temperature these distributions 
start exhibiting two peaks. Notably, the temperature at which the hopping-like motion appears 
for the smaller particles is slightly higher than that at which the hopping-like motion appears for 
the larger ones. We define and calculate a new non-Gaussian parameter whose maximum occurs 
approximately at the time at which the two peaks in the probability distribution of the logarithm 
of displacements are most evident. 

PACS numbers: 61.20 Lc, 64.70 Pf, 61.43 Fs 



I. INTRODUCTION 

Computer simulations have added a great deal to our 
knowledge of relaxation in supercooled liquids. A seminal 
series of studies was performed by Kob and Andersen Q 
who were interested in comparing the results of molecu- 
lar dynamics simulations to the predictions of the mode- 
coupling theory Q . They observed that as the tempera- 
ture was lowered the mode-coupling theory gave a good 
qualitative description of the relaxation of the liquid. 
Specifically, they found that the long-time self-diffusion 
coefficient and a characteristic relaxation time showed 
power law behavior as T c = 0.435 was approached. Since 
their investigation, this temperature has been referred 
to as the mode-coupling transition temperature for the 
Kob- Andersen Lennard- Jones mixture. Later work 
demonstrated that there is no true vanishing of the self- 
diffusion coefficient or true divergence of the a relaxation 
time at T c . This is similar to what is observed experi- 
mentally for most glass formers [{|. There exists power 
law like behavior of the relaxation time as a temperature 
T c is approached and close to T c there is a crossover to a 
different relaxation scenario which extends to the labo- 
ratory glass transition T g , defined as the temperature at 
which the viscosity equals 10 13 poise [MEJ0- 

Many simulations performed since the Kob and An- 
dersen study have resulted in the following relaxation 
scenario in supercooled liquids around the crossover tem- 
perature T c (see Ref. and references therein). The 
motion of particles in a slightly supercooled liquid is ho- 
mogeneous and the self part of the van Hove correlation 
function [8| is approximately Gaussian at all times. At 
lower temperatures the self part of the van Hove corre- 
lation function deviates significantly from Gaussian, the 



motion of the particles is strongly heterogeneous, and, 
on an intermediate time scale (up to so-called exchange 
time), the particles can be separated by their individual 
relaxation times, ft is generally believed that at short 
times the particles are confined to cages of neighbor- 
ing particles and a "cooperative motion" fiTj of parti- 
cles has to occur to facilitate cage escape |l2l Il3j. The 
timescale of this cooperative motion increases with de- 
creasing temperature. In the original studies 0, it 
was assumed that this time scale is around the peak in 
the non-Gaussian parameter ct2(t) — I (<5r 4 ) / (£r 2 ) — 1, 
where Sr is the distance over which a particle moved in 
time t. It has been accepted that the peak position of 
a.2{t) "roughly locates the time of maximum dynamic 
heterogeneity" 0]. While this is a natural first choice, 
there is no a priori reason to choose this time. It has 
been shown that the van Hove correlation function devi- 
ates strongly from a Gaussian distribution over a much 
longer time scale [l^. Moreover, it is also possible to 
identify particles which remain slower than other parti- 
cles for time much longer than the time scale given by the 
peak position of 02 (t) 0. We should point out that, al- 
though the time scale of the peak position of «2 (t) is often 
singled out in connection with dynamic heterogeneity, it 
has also been recognized that there exists dynamic het- 
erogeneity on a much longer time scale, and a four-point 
correlation function that is sensitive to this longer-time 
heterogeneity has been defined and investigated (see Ref. 
[l7j| and references therein). 

In this work we examine the relaxation of a super- 
cooled liquid using Brownian dynamics computer sim- 
ulations. Experiments conducted on Brownian systems 
(i.e. on colloidal suspensions) have so far failed to show 
any significant departure from the power law behavior of 
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the self-diffusion coefficient and the characteristic relax- 
ation time [is|. Thus, it was suggested that the predic- 
tions of the original (also known as "idealized") mode- 
coupling theory provide a good description of the relax- 
ation of the fluid. However, it has since been demon- 
strated using computer simulations that the relaxation 
of a supercooled fluid is the same for Newtonian dy- 
namics m_, stochastic dynamics |4|, and Brownian dy- 
namics |l5| . In particular, in computer simulations de- 
partures from mode-coupling-like power laws and emer- 
gence of hopping-like motion have been observed for all 
three microscopic dynamics. The qualitative difference 
between experimental results (which do not find devia- 
tions from mode-coupling-like power laws) and theoreti- 
cal ones (which do) remains unexplained. 

The focus of this work is the crossover from the high 
temperature homogeneous relaxation of the slightly su- 
percooled fluid to the low temperature relaxation around 
T c . We show that, while the mode-coupling transition 
temperature cannot be unambiguously determined from 
the computer simulation results, the commonly reported 
temperature is the crossover temperature for two differ- 
ent modes of relaxation. Furthermore, we define and 
calculate a time-dependent function that vanishes identi- 
cally for a Gaussian diffusion process. This new function 
7(i) = i ^(5r 2 ) (l/5r 2 ) — 1, hereafter called a new non- 
Gaussian parameter, has a peak at a time at which the 
relaxation seems the most heterogeneous. This time is al- 
ways longer than a characteristic decay time of the inco- 
herent intermediate scattering function [9| that is known 
as the a relaxation time. Moreover, the peak position 
of the new non-Gaussian parameter has temperature de- 
pendence very similar to that of the a relaxation time. 

The paper is organized as follows. In section [n] 
we briefly describe the simulation. In section II I II we 
present results for the mean square displacement, the self- 
diffusion coefficient, the incoherent intermediate scatter- 
ing functions, and the a relaxation time (a preliminary 
report of some of these results appeared in Ref. ^H). We 
show that two different power-law fits can be obtained de- 
pending on the range of temperatures used. In section 
IIVI we present results for the probability distributions 
of the logarithm of single-particle displacements. These 
distributions are sensitive to hopping-like dynamics. In 
section Ivl we define and present the new non-Gaussian 
parameter. In scction lVll we discuss the conclusions that 
can be drawn from this work. 



The simulations are performed with the interaction po- 
tential cut at 2.5 cr Q /3, and the box length of the cubic 
simulation cell is 9.4 uaa- Periodic boundary conditions 
were used. 

We performed Brownian dynamics simulations. The 
equation of motion for the position of the i t h particle of 
type a, rf, is 



= 1f° 
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(1) 



where the friction coefficient of an isolated particle £o = 
1.0 and F" is the force acting on the i t n particle of type 
a, 



3* = -v$TEM 



(2) 



with Vf being the gradient operator with respect to ff. 
In Eq. the dot denotes a time derivative, and the 
random noise ffi satisfies the fluctuation-dissipation the- 
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In Eq. ©, the diffusion coefficient Dq = fcsT/^o where 
kg is Boltzmann's constant and 1 is the unit tensor. 
Since the equation of motion allows for diffusive motion 
of the center of mass, all the results will be presented rel- 
ative the center of mass (i.e. momentary positions of all 
the particles are always relative to the momentary posi- 
tion of the center of mass). We will present the results in 
terms of the reduced units with oaa, ^aa, ^aaI^b, and 
c AA £,a/eAA being the units of length, energy, tempera- 
ture, and time, respectively. The mass of both particles 
are the same and equal to 1.0. 

The equations of motion, Eq. ^ were solved using a 
Heun algorithm with a small time step of 5 x 10~ 5 . We 
simulated the temperatures T = 0.44, 0.45, 0.47, 0.5, 
0.55, 0.6, 0.8, 0.9, 1.0, 2.0, 3.0, and 5.0. We ran a long 
equilibration run (at least half as long as the production 
run) and four production runs at each temperature, ex- 
cept at T = 0.44 where we ran six production runs. The 
results are an average over the production runs, which 
were as long as 6 x 10 s steps long for the lowest temper- 
ature studied. 



II. SIMULATION DETAILS 

We simulated a binary mixture of 800 particles of type 
A and 200 particles of type B that was first considered 
by Kob and Andersen Briefly, the interaction po- 
tential is V a(3 (r) = 4:e a pl(a a p/r) 12 - (cr Q/3 /r) 6 ], where 
a, 13 e {A,B}, e A A = 1-0, o A a = 1-0, <^ab = 1.5, 
&AB = 0.8, €bb = 0.5, and obb — 0.88. The interac- 
tion parameters are chosen to prevent crystallization p| . 



III. MODE-COUPLING-LIKE POWER LAWS 

The mode-coupling theory predicts a power law van- 
ishing of the self-diffusion coefficient and a power law 
divergence of the characteristic relaxation time at the 
mode-coupling transition temperature T c . In many sim- 
ulations and experiments there is a range of temperatures 
where power laws fit the diffusion coefficient and the re- 
laxation time well, and the transition temperature is ob- 
tained from fits of these properties to functions of the 
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FIG. 1: The mean square displacement for the A particles for 
T = 5.0, 3.0, 2.0, 1.5, 1.0, 0.9, 0.8, 0.6, 0.55, 0.50, 0.47, 0.45, 
0.44 listed from left to right. The graph of the mean square 
displacement for the B particles is similar. 
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FIG. 2: Long time diffusion coefficient for the A (closed sym- 
bols) and the B (open symbols) particles. The lines are the 
power law fits to the function a(T — T c ) 7 with T c = 0.435 
fixed. 



form a(T — T c ) 7 . The T c obtained in this manner is gen- 
erally referred to as the mode-coupling temperature. In 
this section, we will present results for the mean square 
displacement, the diffusion coefficient, the self intermedi- 
ate scattering functions and the a relaxation time. More- 
over, we will show that for the system studied in this 
work reasonable power law fits can be obtained for the 
diffusion coefficient and the a relaxation time for a tran- 
sition temperature different from the usually accepted 
mode-coupling transition temperature of T c = 0.435 if a 
different range of temperatures is used for the power law 
fits. 

Shown in Figgis the single-particle mean square dis- 
placement (Sr 2 (t)) = (jr A (t) -r A (0)\ 2X J for the A par- 
ticles. The graph for the B particles looks similar. The 
short time motion is diffusive with a temperature de- 
pendent diffusion coefficient Do = ksT/^o where ks is 
Boltzmann's constant and £o is the friction coefficient 
(recall that in our units ks — 1 and £o = !)• For low 
temperatures a plateau develops in the log-log plot at 
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FIG. 3: Diffusion coefficient for the A (closed symbols) and B 
(open symbols) particles. The solid lines are the "alternative" 

power- law fits to the function a (2) (T - T c (2) ) 7<2) where T c (2) = 
0.401. Dashed lines are fits to the function a(T — T c ) 7 where 
T c = 0.435. 

intermediate times where the mean square displacement 
does not change appreciably. For the lowest temperature 
the plateau region spans several decades in time. 

The plateau region is generally associated with the 
"cage" effect. It represents a localization of the parti- 
cles on intermediate time scales and has been observed in 
many simulations of glassy systems [HQ. Note that there 
is no true plateau in the mean square displacement ver- 
sus time and the slope of the mean square displacement 
versus time decreases monotonically. Therefore there is 
no inflection point of the mean square displacement as 
a function of time. However, there is an inflection point 
in the logarithm of the mean square displacement versus 
the logarithm of time. We use this inflection point to find 
the cage diameter. For T = 0.44 the inflection point oc- 
curs at a value of the mean square displacement of 0.0288 
<jaa for the A particles and 0.0461 <jaa for the B parti- 
cles. These values of the cage diameter corresponds to a 
distance around 0.17 <jaa for the A particles and around 
0.21 o aa for the B particles, which is much smaller than 
the diameter of any particle. After the plateau, the mo- 
tion is again diffusive with a diffusion coefficient D < Dq. 

We determined the long-time self-diffusion coefficient 
D from the slope of (<5r 2 (i)) at long times. The results 
for the A and B particles are shown in figure [21 We ob- 
serve power law behavior similar to what was reported in 
previous simulations of this system using Newtonian 
and stochastic dynamics 0, |j| . Namely, there is power 
law behavior of the diffusion coefficient for temperatures 
between T = 0.8 and T = 0.50. Then there are devia- 
tions from the power law at and below T = 0.47. We fit 
the diffusion coefficients for 0.5 < T < 0.8 for the A and 
B particles to a power law of the form a(T — 0.435) 7 . The 
exponents in the fit are given in the figure, and agree rea- 
sonably well with the exponents found using Newtonian 
and stochastic dynamics 0, 0, 0] . 

We also performed three parameter fits to the diffusion 
coefficient and the a relaxation time for different temper- 
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FIG. 4: Self intermediate scattering functions for the A parti- 
cles for T = 5.0, 3.0, 2.0, 1.5, 1.0, 0.90, 0.80, 0.60, 0.55, 0.50, 
0.47, 0.45, and 0.44 listed from left to right. The correspond- 
ing graph for the B particles is similar. 
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FIG. 5: The a relaxation time for the A (closed symbols) and 
B (open symbols) particles. The lines are fits to the function 
a(T - Tc)' 1 with T c = 0.435 fixed. 



ature ranges |20| . We found that we can obtain good fits 
to D and r a for temperatures 0.44 < T < 0.6. The transi- 
tion temperature depends slightly on the quantity being 
fitted, and ranges from T c = 0.391 to T c = 0.409. We 
found the average transition temperature obtained from 
the fits T c (2) = 0.401 ± 0.009, and then fit each quantity 
to the function a { ^(T - 0.401) 7<2) . We show this final fit 
for the self-diffusion coefficient in Fig. [3J Qualitatively, 
it is clear that the new fit is as good as the standard fit 
shown in Fig. |21 Quantitatively, we evaluate the qual- 
ity of fit by examining X 2 = Ei [{ D ( T i) - Vpi))/*] 2 
where N is the number of data points used in the fit, 
y(T) = a(T — Tc) 1 and Oi is the standard deviation of 
the diffusion coefficient at T^. For each fit we find the 
probability p that the value of x 2 should exceed the calcu- 
lated value by chance given that the model is correct [2l| . 
Higher values of p correspond to a better fit, and the max- 
imum value of p is one. For the transition temperature 
T c (2) = 0.401 and the temperature range 0.44 <T< 0.6, 
p = 0.860 for the A particles and p = 0.902 for the B 
particles. This is significantly better than the values 
p = 0.044 and p = 0.0169 for the A and B particles, re- 
spectively, when the diffusion coefficient was fit over the 
temperature range 0.5 <T< 0.8 and with T c — 0.435. 

We calculated the self intermediate scattering function 
F?(q,t) = ^(Ef=iexp[^-(^(t)-f?(0))]), a G 

{A, B}, shown in Fig.0J for a wave vector q = \q\ around 
the first peak of the partial structure factors for the A 
(q = 7.25) and the B (q = 5.75) particles. The self 
intermediate scattering function decays from its i = 
value of one to zero. For the same temperatures in which 
there is a plateau in the log-log plot of the mean square 
displacement, there is also a plateau in the log-log plot 
of the incoherent intermediate scattering functions. The 
characteristic time for the decay of the self intermediate 
scattering function is the a relaxation time r a , which we 
define as the time when this scattering function is equal 
to 1/e of its initial value. It has been observed that other 
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FIG. 6: The a relaxation time for the A (closed symbols) 
and B (open symbols) particles. The solid line are the "al- 
ternative" power-law fits to the function a^(T — TJ 2 ') ' ' 

(2) 

where Ti J = 0.401. Dashed lines are fits to the functions 
a{T - Tc)'" 1 where T c = 0.435. 



definitions of the a relaxation time results in the same 
temperature dependence. 

The a relaxation time is shown in Fig. |3J Again there 
is power law behavior of the a relaxation time between 
0.5 < T < 0.8, then there are deviations from the power 
law for temperatures at and below T = 0.47. The lines 
in the figure are power law fits to the function a(T — 
0.435) 1 for the temperature range 0.5 < T < 0.8, which 
is the same function and the same temperature range 
as for Fig. [21 The exponents are given in Fig. [S] and 
are close to what was found in simulations of the same 
system using Newtonian 0] and stochastic dynamics 0, 
3 instead of Brownian dynamics. 

As explained above, we also fit the a relaxation time to 
the function a^(T - 0.401)" 7<2> . The result is shown in 
Fig. [5] Again, qualitatively, the new fit is as good as the 
standard fit shown in Fig.[SJ Quantitatively, we evaluated 
the quality of fits using the same procedure as for the 
diffusion coefficient. When T c was set to 0.401 and we fit 
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the a relaxation time for 0.44 < T < 0.6, p = 0.598 for 
the A particles and p — 0.743 for the B particles. When 
the transition temperature was set to T c = 0.435 and we 
used the temperature range 0.5 < T < 0.8, p = 0.15 for 
the A particles and p = 0.35 for the B particles. Thus 
the fits were better for T c (2) = 0.401. 

The values of the scaling exponents obtained from new 
fits for both the self-diffusion coefficient and the a relax- 
ation time are quite a bit greater than those obtained 
from the conventional fit. Also, the difference between 
the exponents for the A and B particles is considerably 
larger. 

We should point out that one could try using a slightly 
different procedure to identify the mode-coupling transi- 
tion temperature. Namely, one could try fit a straight 
line to a plot of, e.g D 1 ^ 1 vs. temperature, where 7 is 
a scaling exponent obtained from solving mode-coupling 
equations. The are two potentials problems with this ap- 
proach: first, it requires solving full time and wave- vector 
dependent mode-coupling equations |22j. Second, using 
this alternative approach one has to neglect the quali- 
tative differences between the predictions of the theory 
and results of the simulations like, e.g. the difference be- 
tween scaling exponents for the A and the B particles and 
the difference between scaling exponents for the diffusion 
coefficients and the relaxation times. 

In conclusion, we find that the mode-coupling temper- 
ature of the Kob- Andersen binary mixture is not unique: 
if a temperature range different from the original one 
is used for fitting the simulation data to the power laws, 
a different mode-coupling transition temperature results. 
In the next section we show that there is a change in the 
dynamics around the commonly reported mode-coupling 
temperature, T c — 0.435. 



IV. PROBABILITY DISTRIBUTIONS OF 
SINGLE PARTICLE DISPLACEMENTS 

Following a procedure suggested previously 
we examined the probability distributions of the loga- 
rithm of single-particle displacements P(\og 10 (Sr); t) at a 
time t. Multiple peaks in P(log 10 (Sr); t) at a time t' have 
been clearly observed in simulations of model colloidal 
gels Q , and provide evidence of populations of fast and 
slow particles at t'. Note that P(log 10 (Sr);t) is defined 
such that the integral / 1 P(x; t)dx is the fraction of par- 
ticles whose value of log 10 (i5r) is between xo and x\. Fur- 
thermore, the probability distribution P(log 10 (Sr); i) can 
be obtained from the self van Hove correlation function 
@ since P(log 10 (c5r); t) = ln(10)47rfr 3 G s (5r, t). There 
are some general properties of P(\og w (Sr); t) if the van 
Hove correlation function is Gaussian. If the motion 
of a tagged particle is diffusive at all times with a dif- 
fusion coefficient D, then the self van Hove correlation 
function G s (Sr,t) = (l/(47r£»i)§) exp(-Sr 2 /ADt). For a 
Gaussian van Hove function, the shape of P(log 10 (<5r); t) 
is independent of time [25j, and the peak height of 



P(log w {Sr);t) = log e (10)v / 54/V e 3 / 2 » 2.13. Devia- 
tions from this height represents deviations from Gaus- 
sian behavior of G s (r, t) . 

Shown in Fig. is P(\og 10 (5r); t) for the A particles 
at T = 5.0, 0.6, 0.47, and 0.44. For each temperature 
several different times are shown. The thicker lines corre- 
spond to later times. At T = 5.0 there is little deviation 
from Gaussian behavior: there is only one peak whose 
height is close to 2.13 at all times. The peak position 
moves to larger distances for larger times. For T = 0.6, 
there are deviations from Gaussian behavior manifested 
in the reduced height of the peak of the distribution, but 
there is still only one peak for all times. At T = 0.47, 
the deviations from Gaussian behavior are stronger and 
the distribution becomes very broad at a time which cor- 
responds to right after the plateau of the mean square 
displacement. At T = 0.44 there are two distinct peaks. 
The position of the second peak depends on time, but 
when the height of both peaks are approximately equal 
the position of the second peak is around log 10 (<5r) = 0, 
thus Sr = 1.0 = o~aa- 

The probability distributions P(\og 10 (5r);t) provide 
clear evidence that there are populations of particles with 
different mobilities. These probability distributions are 
similar to the distributions observed in model colloidal 
gels close to a gelation transition [23| and in dense sys- 
tems with purely repulsive interactions The shape 
of the distributions strongly suggests a heterogeneous 
hopping-like motion for at least a fraction of the par- 
ticles. We find that the hopping rate of the particles 
varies greatly between particles, and the typical size of 
the hopping length for the A (larger) particles is equal to 
one particle diameter. 

The two peaks in P(\og 10 (5r); t) are clearly defined at 
a higher temperature for the B particles, Fig. El than for 
the A particles. For T = 5.0, P(\og 10 (Sr); t) for the B 
particles is similar to what is observed for the A particles. 
The distribution P(\og 10 (Sr); t) broadens and the peak 
height decreases significantly for the B particles start- 
ing at T = 0.6, and we observe the two peak structure 
starting at T = 0.5. The two peaks are well defined for 
T = 0.47 and are very prominent for T = 0.44. At these 
low temperatures there is a clear distinction between mo- 
bile and immobile particles. The second peak occurs 
around log 10 (£r) « 0.086 for the B particles, which cor- 
responds to a value of Sr « 1.25. This suggests that the 
typical length of the particle jumps are slightly larger 
for the B particles than for the A particles. Note that 
there are small secondary peaks in P(log 10 (Sr); t) for the 
B particles for the longer times. This is expected if the 
activated hopping is the dominant relaxation process. 

The size dependence of the distributions of the dis- 
placements of the particles have been observed in simu- 
lations of colloidal gels by Puertas et al. [23| . They simu- 
lated a system of soft core polydisperse particles with an 
average particle radius a. They noticed that the distri- 
butions of squared displacements were more bimodal for 
particles with a smaller radius (see Fig. 6 in Ref. |23|). It 
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FIG. 7: The probability of the logarithm of single particle displacements P(log 10 (Sr);t) for the A particles. The wider lines 
indicate increased time t. a: T = 5.0 (a) for t = 0.01, 0.05, 0.1, 5 and 10 listed from left to right, b: T — 0.60 for t = 1.0, 5, 15, 
25 and 30 listed from left to right, c: T = 0.47 for t = 100, 200, 300, 400, 500, and 800 listed from left to right, d: T = 0.44 
for t = 1000, 2000, 3000, 4000, 5000, 6000, and 10000 listed from left to right. 



is likely that the size dependence of dynamic heterogene- 
ity is a general feature of systems with slow dynamics. 

In earlier studies [14 . |23( probability distributions of 
particle displacements were examined at a time t* such 
that (oY 2 (i*)) = 10a 2 , where a is a measure of the av- 
erage size of the polydisperse particles. While bimodal 
distributions were observed for this t* , there was no clear 
justification for selecting this particular timescale. In the 
next section we introduce a new non-Gaussian parame- 
ter whose peak position allows us to identify the time 
in which the two peaks in P(\og w (8r);t) are around the 
same height. Moreover, we will demonstrate that the 
peak position of the commonly used non-Gaussian pa- 
rameter Q!2(i) does a poor job of identifying this time. 



V. NEW NON-GAUSSIAN PARAMETER 

Many simulations which examine heterogeneous dy- 
namics focus on a timescale which is found from 
the peak of the non-Gaussian parameter oli (t) = 

| (Jr 4 (£)) / (Sr 2 (t)) 2 — 1. It has been observed that close 
to T c = 0.435, a 2 (i) has a peak occurring at the so- 
called late (3 regime, i.e. some time before the a relax- 
ation time. Moreover, for T < 0.8 with decreasing tem- 



perature the value of the non-Gaussian parameter at the 
a relaxation time, 0:2 (tv>), is a decreasing fraction of its 
maximum value |l5j . In particular, at T = 0.44 the non- 
Gaussian parameter at the a relaxation time is approxi- 
mately equal to 0.62 and 0.35 of its maximum value for 
the A and B particles, respectively. We observe that the 
two peaks in P(log 10 (<5r);t) (for temperatures in which 
the two peaks are clearly defined) are about the same 
height at some time after the a relaxation time. Thus 
ai (t) is small compared to its maximum value at a time 
when the two peaks are very prominent. 

We propose to use a new function that, as does a 2 (t), 
quantifies deviations from a Gaussian distribution of dis- 
placements, 

7(*) = (^) ~ 1= (4) 

where Sr is the distance over which the particle moved 
in time t. We will show that this new non-Gaussian pa- 
rameter identifies the time in which the two peaks in the 
probability distribution of the logarithm of single-particle 
displacements are most evident. The parameter 0:2 (i) is 
significantly influenced mostly by particles which move 
farther than expected from a Gaussian distribution of 
particle displacements. In contrast, the parameter j(t) 
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FIG. 8: The probability of the logarithm of single particle displacements P(log 10 (Sr);t) for the B particles. Wider lines 
indicates increased time t. (a): T = 0.5.0 for t = 0.01, 0.05, 0.1, and 0.5 listed from left to right, (b): T = 0.60 for t = 1, 5, 10, 
15, 20 and 25 listed from left to right, (c): T = 0.47 for t = 100, 200, 300, 400, and 500 listed from left to right, (d): T = 0.44 
for t = 600, 700, 800, 1000, 3000, 4000, and 5000 listed from left to right. 



weights very strongly the particles which have not moved 
as far as expected from a Gaussian distribution of parti- 
cles displacements through the term (l/<5r 2 ) and weights 
the particles which move farther than expected from a 
Gaussian distribution through the (<5r 2 ) term. The fac- 
tor of 1/3 ensures that 7(i) is zero when the self part of 
the van Hove correlation function is Gaussian. 

Shown in Figgis j(t) for the A and B particles. At 
short times, the motion of the particles are Gaussian and 
7(i) is close to zero. At intermediate times there is a 
peak in "f(t) whose height and position increases with 
decreasing temperature. The peak in j(t) is larger for 
the B particles, which is what we expect from examining 
the self part of the van Hove correlation function and 
P(log 10 (<5r); t). At long times j(t) decays to zero. 

It is interesting to note the temperature dependence 
of the the average mean square displacement at the peak 
position of j(t): for T > 0.8 the average mean square dis- 
placement is independent of temperature and it is around 
0.23(T AA for the A particles and 0.31a AA for the B par- 
ticles. For T < 0.8, the average mean squared displace- 
ment at the peak position of j(t) increases with decreas- 
ing temperature and at T — 0.44 it reaches the values of 
approximately 0.84a AA for the A particles and 2A0a AA 
for the B particles. It should be pointed out that the 



temperature at which the average mean squared displace- 
ment at the peak position of j(t) starts increasing is close 
to the so-called onset temperature identified for the Kob- 
Andersen model by Brumer and Reichman |27j . 

The insert to each graph in Fig.^is the peak position of 
the new non-Gaussian parameter, j(t), denoted by T nng 
compared to the a relaxation time and the peak position 
of the commonly used non-Gaussian parameter, a%{€), 
denoted by T ng . Notice that the peak position of 7(f) is 
always greater than the a relaxation time, but it has the 
same temperature dependence. This is in contrast to r ng 
which is greater than the a relaxation time and equal to 
Tnng at higher temperatures, but increases slower with 
decreasing temperature than r Q and T nng . For the lowest 
temperatures, r ng is much smaller than either r Q and 

Tnng • 

The peak position approximately corresponds to the 
time in which the two peaks of P(log 10 (5r) ; t) are of equal 
height. Shown in Fig.HUIand Fig.HTlis the P(log 10 (<5r); t) 
for the A and B particles at the peak position of a2(t) 
(dashed lines) and j(t) (solid lines) for temperatures of 
T = 0.47 and T = 0.44. Notice that at the time corre- 
sponding to the peak position of ui (t), there is at most 
a shoulder in P(\og 10 (Sr); t). For the time which corre- 
sponds to the position of the peak in j(t), the two peaks 
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FIG. 9: The new non-Gaussian parameter ~/(t) for T — 5.0, 
3.0, 2.0, 1.5, 1.0, 0.9, 0.8, 0.6, 0.55, 0.50, 0.47, 0.45, 0.44 for 
the A (a) and B (b) particles. The insets are the position 
of the peak of j(t), (open circles) compared to a relaxation 
time (closed circles) and the peak position of ai (t) (squares) 
versus temperature. 

in P(log 10(<5r); t) are of similar height and the definition 
of mobile and immobile particles is clear. For compari- 
son, the self part of the van Hove correlation function at 
the peak position of (X2(t) (dashed lines) and j(t) (solid 
lines) are shown as inserts in Fig. 1 101 and Fig. 1111 

VI. CONCLUSIONS 

The mode-coupling theory has been used extensively 
to describe the relaxation in supercooled liquids and the 
glass transition. It correctly describes many qualitative 
features of the glass transitions observed experimentally 
and in computer simulations. The most notable success 
of the mode-coupling theory is that it correctly describes 
the two step decay of the intermediate scattering func- 
tions and the qualitative features of the mean square dis- 
placement as a function of time. 

The idealized mode-coupling theory predicts power law 
divergence of the a relaxation time and power law van- 
ishing of the self-diffusion coefficient at a transition tem- 
perature T c . It has been observed that there is power 
law like behavior of the a relaxation time and the diffu- 
sion coefficient in simulations and experiments close to a 



FIG. 10: The probability of the logarithm of the displace- 
ments at the the peak position of the the old non-Gaussian 
parameter a.2{t) (dashed line) and the new non-Gaussian pa- 
rameter 7(t) (solid line) at T = 0.47 for the A particles (a) 
and the B particles (b). The inserts are the self part of the 
van Hove correlation functions at the peak position of a2(t) 
(dashed line) and y(t) (solid line). 

temperature T c . However, the structural arrest predicted 
by the mode-coupling theory does not occur at T c , but 
rather there appears to be a crossover to a different relax- 
ation scenario and the mode-coupling transition is said 
to be "avoided" . 

We performed Brownian dynamics simulations of a fre- 
quently studied glass forming binary mixture for a large 
range of temperatures. The temperature dependence of 
the a relaxation time and the diffusion coefficient were 
similar to what was observed in previous simulations of 
the same system 0, 0, Q . We found that the diffusion 
coefficient and the a relaxation time can also be fit to 

(2) 

a power law where the transition temperature is T c = 
0.401 which is lower than the generally accepted mode- 
coupling transition temperature T c — 0.435. However, 
by examining the probability distributions of the loga- 
rithm of single particle displacements, P(\og w (Sr); t), we 
demonstrated that the crossover from the high temper- 
ature diffusive relaxation of the particles to low temper- 
ature hopping-likc motion occurs near the generally ac- 
cepted mode-coupling transition temperature T c — 0.435. 

The change in the relaxation processes are evident 
when one examines the distribution of the logarithm of 
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FIG. 11: The probability of the logarithm of the displace- 
ments at the peak position of the old non-Gaussian parame- 
ter 02 (t) (dashed line) and the new non-Gaussian parameter 
y(t) (solid line) for T — 0.44 for the A particles (a) and the 
B particles (b). The inserts are the self part of the van Hove 
correlation functions at the peak position of 02 (i) (dashed 
line) and y(t) (solid line). 



single particle displacements P(\og w (Sr); t). At higher 
temperatures there is one peak in P (\og 10 (5 r);t) whose 
position increases for increased time. For lower temper- 
atures, the distribution becomes broad at a time scale 
right after the plateau region of the mean square dis- 
placement. At the lowest temperatures examined in this 
study, there are two distinct peaks in P(\og 10 (Sr); t) for 
both the A and B particles. The two peaks are evidence 
that on an intermediate time scale the particles can be 



separated by their individual relaxation time, and thus 
are dynamically heterogeneous. Moreover, the minimum 
between the peaks is smaller for the B particles than for 
the A particles at a fixed temperature. The dependence 
on particles size of the probability distributions is similar 
to what was observed in earlier simulations [Til |23| . 

One possible interpretation of our results could be that 
the mode-coupling temperature does not have any phys- 
ical significance and the crossover in supercooled liquid's 
dynamics is very smooth |29| . While such an interpre- 
tation cannot be excluded, we would like to advocate a 
more cautious conclusion: in order to identify a crossover 
temperature (or a narrow crossover temperature range) 
that one could interpret as mode-coupling temperature, 
one has to investigate not only power-law fits to trans- 
port coefficients and/or relaxation times but also whether 
microscopic dynamics is homogeneous and diffusive-like 
or heterogeneous and hopping-like. We would like to em- 
phasize that although the mode-coupling theory, in its 
standard form, cannot describe hopping-like dynamics, 
it does provide a reasonable |22j description of dynamics 
in moderately supercooled fluids. 

We found that the typical non-Gaussian parameter 
ct2(t) does a poor job of identifying the time scale on 
which the heterogeneous, hopping-like motion is most 
evident. We defined a new non-Gaussian parameter 
7(i) = I (<5r 2 ) (l/<5r 2 ) — 1. For temperatures in which 
there are two peaks in P(\og 1Q (6r);t), ^(t) has a peak 
occurring at a time T nng in which the two peaks have 
approximately the same height. For times in which only 
one peak is present, the peak position of j(t) identifies 
the time in which P(log 10 (<5r); t) is widest. The position 
of the peak for at2(t) and j(t) are the same for high tem- 
peratures, but the peak position of 0:2 (i) increases slower 
with decreasing temperature than the peak position of 
j(t). Moreover, the peak position of j(t) has the same 
temperature dependence as the a relaxation time. 
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